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affects  the  DLD  of  the  langasite  resonator.  The  DLD  measurement  results  for  langasite  have  been  compared  with  literature 
values  for  quartz,  langasite,  and  langatate,  and  with  additional  new  measurements  for  a  GaP04  resonator  of  type  R-30  (-11.1° 
rotated  Y-cut).  We  have  also,  for  the  first  time,  performed  an  uncertainty  analysis  for  the  measured  drive  level  sensitivity.  The 
result  has  led  to  a  journal  paper  in  [9]. 

3.  Analysis  of  a  surface  acoustic  wave  quartz/langasite  resonator 

Stress-induced  surface  acoustic  wave  velocity  shifts  are  analyzed  for  langasite  (LGS)  SAW  resonators.  The  analytical 
methodology  has  been  verified  by  comparing  experimental  results  and  analytical  results  for  quartz  resonators.  LGS  SAW 
resonators  with  Euler  angles,  which  are  most  sensitive  and  least  sensitive  to  diametrical  forces,  are  determined  and  their 
applications  in  force  sensors  and  resonators  with  minimum  acceleration  sensitivity  are  discussed.  Sensitivity  of  the  analytical 
results  to  different  groups  of  nonlinear  material  constants  is  discussed;  it  was  found  that  for  specific  configurations,  failure  to 
include  the  third-order  piezoelectric  constants,  dielectric  constants  and  electrostrictive  constants  may  lead  to  a  significant 
calculation  error.  Surface  acoustic  waves  propagating  on  an  LGS  square  plate  subject  to  bending  moment  along  the 
propagation  direction  and  normal  to  the  propagation  direction  are  analyzed;  it  was  found  that  the  average  moment-induced 
velocity  shift  of  LGS  resonators  are  comparable  to  quartz  resonators.  Analyses  of  the  sensitivity  of  the  results  to  different 
groups  of  nonlinear  material  constants  show  that  for  some  specific  wave  propagation  directions,  failure  to  include  the  third-order 
piezoelectric  constants,  dielectric  constants,  and  electrostrictive  constants  may  lead  to  large  errors.  The  result  has  led  to  a 
journal  publication  in  [10]. 

4.  Analysis  of  multi-layer  Aluminum  Nitride  thin  film  resonators 

We  have  performed  a  theoretical  analysis  for  a  two-layered  and  mulitlayer  zig-zag  c-axis  inclined  zinc  oxide  thin  film  resonator. 
We  derived  the  frequency  equation  and  mode  shape  for  these  resonators  based  on  the  linear  piezoelectric  theory.  The 
impedance  characteristics  of  the  Thin  Film  Bulk  Acoustic  Wave  Resonator  (FBAR)  are  derived  and  compared  with  previous 
experimental  results.  The  result  has  led  to  two  journal  papers  in  [11]  and  [12].  The  application  of  our  analysis  led  to  another 
journal  paper  on  automotive  oil  viscosity  sensor  [13]. 

5.  Analysis  of  wide-band  piezoelectric  energy  harvesters 

T raditional  piezoelectric  energy  harvesters  operate  in  a  single  resonant  frequency.  If  the  excitation  frequency  is  out  of  the 
operation  frequency  range,  the  energy  harvester  will  fail  to  generate  meaningful  energy.  We  have  developed  a  new  multi-beam 
piezoelectric  energy  harvester  to  overcome  this  barrier.  The  energy  harvester  we  proposed  features  8  times  wider  operational 
frequency  bandwidth  and  has  an  immediate  impact  on  piezoelectric  energy  harvester  design.  The  result  led  to  two  journal 
publications  in  [14]  and  [15]. 

6.  Free  vibration  of  quartz  resonators 

New  results  for  a  better  understanding  of  the  design  of  a  thickness-shear  mode  AT-cut  quartz  crystal  resonator  were  obtained 
by  studying  a  series  of  free  vibration  problems. 


The  specific  problems  analyzed  are  summarized  below. 


A  rectangular  trapped  energy  resonator  partially  electroded  at  the  center  with  rectangular  electrodes  was  analyzed  in  [16].  A 
single-mode,  three-dimensional  equation  governing  the  thickness-shear  displacement  was  used.  A  Fourier  series  solution  was 
obtained.  The  effects  of  the  electrode  size  and  thickness  on  the  trapped  modes  were  examined.  Results  show  the  vibration  of 
the  trapped  thickness-shear  modes  is  mainly  under  the  electrodes  and  decays  rapidly  outside  the  electrodes. 

Shear-horizontal  free  vibrations  of  an  elastic  cylinder  with  an  oblate  elliptical  cross  section  and  a  traction-free  surface  were 
studied  in  [17].  It  can  represent  a  contoured  resonator.  The  exact  vibration  modes  and  frequencies  were  obtained.  The  energy¬ 
trapping  behavior  of  these  modes  was  examined,  and  the  results  show  trapped  thickness-shear  and  thickness-twist  modes. 

Straight-crested  waves  and  vibration  modes  with  variations  along  the  x3  direction  only  in  an  AT-cut  quartz  plate  resonator  near 
the  operating  frequency  of  the  fundamental  thickness-shear  mode  were  analyzed  in  [18].  Previous  researchers  mainly  analyzed 
xl-dependent  modes.  Mindlin’s  two-dimensional  equations  for  anisotropic  crystal  plates  were  used.  Dispersion  relations  and 
frequency  spectra  of  the  five  relevant  waves  were  obtained.  Certain  known  values  of  the  plate  length/thickness  ratio  needed  to 
be  excluded  to  avoid  unwanted  couplings  between  the  resonator  operating  mode  and  other  undesirable  modes.  Research 
showed  that  an  additional  series  of  discrete  values  of  the  plate  length/thickness  ratio  also  needed  to  be  excluded.  For  overtone 
modes,  the  exact  frequency  spectra  for  x3-dependent  modes  were  obtained  in  [19]. 

A  series  of  new  theoretical  results  on  thickness-shear  and  thickness-twist  frequencies  and  modes  were  obtained  using  the 
scalar  differential  equations  by  Tiersten-Smythe  and  Stevens-Tiersten.  These  results  include  the  solutions  of  a  circular 
resonator  [20],  a  rectangular  trapped  energy  resonator  [21],  and  two  different  contoured  resonators  [22],  [23]. 

7.  Electrically  forced  vibration 

A  theoretical  analysis  on  electrically  forced  vibration  of  the  thickness-shear  vibration  of  a  resonator  with  couplings  to  both 
flexure  and  face-shear  modes  was  performed  in  [24].  Mindlin’s  two-dimensional  equations  for  piezoelectric  plates  were 
employed.  Electrically  forced  vibration  solutions  were  obtained  for  three  cases:  pure  thickness-shear  mode  alone;  two  coupled 
modes  of  thickness  shear  and  flexure;  and  three  coupled  modes  of  thickness  shear,  flexure,  and  face  shear.  The  admittance 
was  calculated  and  its  dependence  on  the  driving  frequency  and  the  length/thickness  ratio  of  the  resonator  was  examined. 
Results  show  that  near  the  thickness-shear  resonance,  the  admittance  assumes  maxima,  and  that  for  certain  values  of  the 
length/thickness  ratio,  the  coupling  to  flexure  causes  severe  admittance  drops,  while  the  coupling  to  the  face-shear  mode 
causes  additional  admittance  changes  previously  unknown  and  unconsidered  in  current  resonator  design. 

8.  Transient  effects 

Most  of  the  theoretical  results  on  quartz  resonators  in  the  literature  are  for  time-harmonic  vibrations.  We  analyzed  electrically 
forced,  transient  thickness-shear  vibration  of  a  quartz  resonator  [25].  A  trigonometric  series  solution  was  obtained  with  which  a 
few  transient  effects  in  a  piezoelectric  resonator  under  a  sudden  change  or  fluctuation  of  the  driving  voltage  amplitude  or 
frequency  were  calculated  and  examined.  It  was  found  that,  for  a  1 .65  MHz  resonator  with  a  material  Q=105,  under  a  sudden 
change  of  the  driving  voltage  amplitude,  the  transition  of  the  resonator  vibration  takes  about  0.1  s,  and  that  a  10  Hz  deviation  of 
the  driving  frequency  from  the  fundamental  resonance  causes  a  significant  drop  of  the  vibration  amplitude.  The  analysis  in  [25] 
was  complicated.  A  much  simpler  method,  which  was  developed  later  in  [26],  can  describe  the  amplitude  evolution  of  transient 
thickness-shear  vibrations  of  a  quartz  resonator  efficiently  and  accurately. 

9.  Second-order  in-plane  acceleration  sensitivity 

We  also  performed  a  theoretical  analysis  on  the  second-order  in-plane  acceleration  sensitivity  of  a  Y-cut  quartz  thickness-shear 
mode  resonator.  The  second-order  nonlinear  theory  of  elasticity  for  anisotropic  crystals  was  used  to  determine  the  biasing  fields 
in  the  resonator  under  in-plane  acceleration.  The  acceleration-induced  frequency  shift  was  determined  from  a  perturbation 
analysis  based  on  the  plate  equations  for  small-amplitude  vibrations  superposed  on  a  finite  bias.  For  a  structurally  symmetric 
resonator  under  in-plane  acceleration,  although  its  first-order  acceleration-induced  frequency  shift  is  known  to  be  zero,  its 
second-order  frequency  shift  was  found  to  be  nonzero  and  is  quadratic  in  the  acceleration.  An  estimate  of  the  second-order 
frequency  shift  is  found  to  be  of  the  order  of  =10-18,  10-16,  and  10-14  when  the  acceleration  is  1  g,  10  g,  and  100  g, 
respectively.  The  result  will  be  submitted  to  a  journal  for  publication. 

10.  The  developed  software  package 

The  “Crystal  resonator  design  and  analysis  system”  software  package  is  based  on  MATLAB  programming  language.  The 
software  requires  MATLAB  201 1  or  above  to  run.  This  software  package  helped  shorten  the  design  cycle  of  a  military  electronic 
system  by  improving  the  resonator  design  efficiency.  The  software  has  three  major  modulus-Material  transform  modulus, 
“frequency  and  mode  shape  analysis  modulus”  and  “frequency  stability  analysis  modulus”. 


1 1 .  Assessment  of  langasite  for  electronic  system  components 


We  have  measured  extensively  some  nonlinear  linear  effect  for  langasite  crystal  resonators  mentioned  in  Section  2  of  this 
report.  These  nonlinear  effect  includes  force  frequency  effects,  electroelastic  effect  and  drive  level  dependence.  Here  we 
provide  an  assessment  of  this  material  for  electronic  system  component  based  on  our  measurement  results: 

a.  For  the  force-frequency  effect  measurements,  we  have  concluded  that  the  langasite  resonators  have  a  relatively  lower  force- 
frequency  effect  compared  with  quartz  resonators,  which  makes  them  a  better  candidate  for  stress  compensated  resonators 
which  has  lower  phase  noise.  This  type  of  resonators  is  ideal  for  vehicle-borne  or  airborne  electronic  device  such  as  GMTI 
(Ground  Moving  Target  Indicator).  The  Minimum  Detectable  Velocity  (MDV)  of  GMTI  is  limited  by  phase  noise  of  the  crystal 
oscillators  mounted  on  it.  With  the  langasite  resonator,  the  MDV-one  of  the  major  performance  benchmark  of  a  GMTI,  could  be 
improved  significantly.  In  addition,  langasite  resonators  could  extend  the  limit  of  airborn  communication  systems.  Since  the 
noise  from  a  transmitter  in  one  channel  extends  to  neighboring  channels,  as  the  number  of  transmitter  grows,  the  noise  can 
accumulate  to  the  point  at  which  receivers  can  no  longer  function  properly.  For  instance,  the  noise  from  a  typical  commercial 
quartz  oscillator  (with  2x  /g  vibration  sensitivity)  limits  the  number  of  users  to  less  than  100  per  transponder.  Low-phase  noise 
langasite  resonators  can  allow  as  many  as  1200  users  per  transponder.  However,  although  we  have  seen  the  promising 
application  of  langasite  resonators  to  military  and  civil  electronic  systems,  we  have  not  been  able  to  measure  the  acceleration 
sensitivity  and  characterize  the  phase  noise  yet,  and  we  have  not  measured  the  influence  of  vibration  to  the  resonant 
frequency.  If  further  funding  is  available,  we  will  measure  the  acceleration  sensitivity  and  phase  noise  of  langasite  resonators, 
which  will  make  us  move  significantly  toward  the  real  application  in  electronic  systems.  Langasite  also  has  a  very  good  high 
temperature  performance  as  its  piezoelectricity  will  remain  till  its  melting  point  around  1400  deg  C,  which  make  it  very  promising 
for  high  temperature  sensing  application  such  as  pressure  and  temperature  sensing  in  harsh  environments.  The  application 
aspects  are  broad  such  as  the  monitoring  of  jet  engines  and  aerospace  propulsion  system,  or  the  monitoring  of  coal/gas-fired 
power  plants,  et.  Al. 

b.  For  the  electroelastic  effect  (electric  field-frequency)  measurements,  the  most  important  discovery  is  that  we  find  that  the 
langasite  resonators  generally  have  one-order-of -magnitude  larger  electric  field-frequency  effect  compared  with  traditional 
quartz  resonators.  This  new  discovery  will  enable  a  highly  sensitive  MEMS  resonant  electric  field  sensor  that  has  the  benefits  of 
high  stability,  low  cost,  low  hysteresis,  low  power  and  a  simple  structure.  It  can  be  used  for  the  measurement  of  both  DC  and 
AC  electric  fields  with  high  sensitivity. 

c.  For  the  drive  level  dependence  (DLD)  measurements,  the  measured  langasite  resonators  exhibit  higher  DLD  compared  with 
quartz  resonators,  further  experiments  and  modelling  for  DLD  of  langasite  resonators  with  other  cut  angles  are  needed  to  make 
a  more  solid  conclusion. 
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Technology  Transfer 


Introduction  of  the  Crystal  resonator  design  and  analysis  system  software 


The  “Crystal  resonator  design  and  analysis  system”  is  software  based  on  MATLAB 
programming  language.  The  software  requires  MATLAB  2011  or  above  to  run.  This  software 
package  helped  shorten  the  design  cycle  of  a  military  electronic  system  by  improving  the 
resonator  design  efficiency.  The  software  has  three  major  modulus-Material  transfonn  modulus, 
“frequency  and  mode  shape  analysis  modulus”  and  “frequency  stability  analysis  modulus”.  The 
material  transfonn  modulus  converts  the  material  constants  including  second-order  and  third- 
order  elastic  constants,  piezoelectric  constants,  dielectric  constants,  and  electrostrictive  constants 
to  a  new  coordinate  system.  By  simply  inputting  the  rotation  angle  <t>,  0,  T,  user  would  be  able  to 
calculate  the  new  second-order  and  third-order  elastic  constants  for  both  quartz  and  langasite 
single  crystals,  the  detailed  introduction  is  in  a  separate  file-Instruction-Material  constant 
transformation.docx. 

The  frequency  and  mode  shape  analysis  modulus  calculate  the  resonant  frequency  and  plot  the 
mode  shape  for  both  bulk  acoustic  wave  (BAW)  and  surface  acoustic  wave  (SAW)  resonators.  It 
included  three  sub-modulus-“SC  cut  analysis  modulus”,  “AT  cut  analysis  modulus”  and 
“Rayleigh  wave  analyses  modulus”.  The  SC  cut  analysis  modulus  allow  users  to  define  the 
resonator  geometry  for  either  rectangular  or  circular  shape,  it  also  provided  electrode  or 
unelectroded  options.  After  defining  resonator  geometry  and  mode  order,  the  user  could  plot  the 
3-D  resonator  mode  shape  and  resonant  frequency.  The  detailed  introduction  can  be  found  in  a 
separate  file-Instruction-Contourresonator.docx. 

Similar  to  the  SC  cut  analysis  modulus,  the  AT  cut  analysis  modulus  allow  users  to  define 
resonator  shape,  geometry,  electrode  condition  and  plot  resonator  mode  shape  and  detennine  the 
resonant  frequency,  the  detailed  introduction  can  be  found  in  a  separate  file-Instruction- 
Contourresonator.docx. 

The  Rayleigh  wave  analysis  modulus  performs  the  saw  wave  frequency  or  wave  length 
calculations  and  plot  the  related  wave  mode.  It  allows  user  to  select  wave  frequency  calculation 
or  wave  length  calculations.  The  detailed  theoretical  background  can  be  found  in  a  separate  file- 
Instruction-Rayleighwave.docx 

The  frequency  stability  analysis  modulus  includes  three  sub-modulus-acceleration  sensitivity 
modulus,  temperature  effect  modulus  and  electric  field  frequency  modulus.  The  acceleration 
sensitivity  modulus  allows  users  to  define  the  resonator  geometry,  mode  order,  and  material 
orientation  and  calculate  the  acceleration  sensitivity,  the  detailed  theoretical  background  can  be 
found  in  a  separate  file-  Instruction-accelerationsensitivity.docx. 

The  temperature  frequency  and  electric  field  frequency  modulus  will  be  developed  in  future. 


Software  installation 


Firstly,  set  up  a  new  folder  in  C  drive,  the  name  of  the  folder  should  be  CRDV 1  . 
C:\CRDV1 

Secondly,  please  copy  all  fdes  into  the  folder  (C:\CRDV1) 

Thirdly,  click  the  Controlplatform.fig  to  run  the  whole  program. 


Material  constant  transformation 


The  foundation  for  the  material  constants  transformation  is  based  on  tensor  transformation 
principle. 


tkmp..  QlkQnmQqp  •••  tlnq... 

While  Q  is  the  material  rotation  matrix.  The  second-order  elastic,  piezoelectric  and  dielectric 
constants  were  transformed  to  new  coordinate  (Fig.  1 )  and  are  shown  in  the  front  panel.  The 
third-order  elastic,  piezoelectric,  dielectric  constants  and  electrostricitive  constants  are 
transfonned  to  the  new  coordinate  as  well.  Due  the  large  number  of  these  nonlinear  constants, 
the  transformed  constants  are  stored  in  a  file  “mydata.mat”  when  click  the  “save”  button  in  the 
front  panel.  Users  can  access  all  transformed  constants  by  loading  this  file  in  the  MATLAB 
command  window. 


Fig.  1.  Plate  orientations. 


Instructions 


2  a 

Figure  1 .  A  simply  supported  rectangular  plate  under  normal  acceleration. 

Purpose:  Calculating  the  displacement  gradient  field  of  the  biasing  deformation  in  a  simply 
supported,  AT-cut  quartz  plate  under  normal  acceleration. 

Input: 

a  — half  of  the  plate  length  (unit:  mm), 
b  —  half  of  the  plate  width  (unit:  nun), 
h  —  half  of  the  plate  thickness  (unit:  mm), 
ai  —  normal  acceleration  (unit:  m/s2). 

Output: 

Wij —  displacement  gradient  (unit:  10"6),  i,j  =1,  2,  3. 


The  deflection  w^ix ,,jc3)  of  a  simply  supported  plate  induced  by  normal  acceleration  is 
governed  by  [1]: 

—h  7 ABCD^ifiDAB  2pha2  =  0 ,  A,  B,  C,  D  — 1,  3.  (1) 

The  boundary  conditions  are: 

wf  =  0  ,  K\\ }  =  0  ,  at  Xj  =  ±a ,  (2) 

wf  =  0  ,  =  0 ,  at  x3  =  ±b .  (3) 

The  solution  can  be  expressed  by  Fourier  series  [1]: 


where 


w2]  =  I]  I]4»„cos(amx1)cos(K-„x3), 

m= 1,3,5  n= 1,3,5 


m n  yitc  .  .  _ 

,  *■„=  — ,  m,  n= 1,3,5,..., 
la  lb 
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mna2h2[yna4m+y33K:*  +(2^13  +4^55)«2a:„4] 


(6) 


We  note  that  the  Amn  in  [1]  has  an  error  and  (6)  is  believed  to  be  correct.  Once  w{20)  is 
known,  we  have  [1] 


Thus  the  plate  strains  E 


Other  plate  strains  Ejp 
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cn  are  given  by  [1]: 
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are  given  by  [1]: 


(7) 


(8) 


E"}  =  -4vcvsEf  ,R,S=\,3,5,W,V=2,4,6,  (9) 

where  R,  S,  W,  V  are  indices  in  compact  index  notation.  The  three-dimensional  strain  field  can 
be  obtained  based  on  (8)  and  (9)  [1]: 


E  =x  E{1) 

C'KL  a2^KL  ’ 

where  K,L=  1,  2,  3.  The  anti-symmetric  part  of  the  displacement  gradient  is: 

According  to  [1],  QA7  can  be  further  simplified  as: 


(10) 


(11) 


D.2A  —  w2  J  +  X2E2A  ’  ^13  _  0  . 


(12) 


Finally  the  displacement  gradients  are  given  by: 

WK,L  =EKL+^LK- 

For  convenience,  we  list  the  displacement  gradient  of  an  AT-cut  plate  as  follows: 

00  00 

I  2X  «3cos(«„,x1)cos(x-„x3), 


(13) 


Wu  =  *2 


m= 1,3,5  n= 1,3,5 


Wi.2  =°.l' 744^2  I]4,m«',AAin(a,,Ai)sin(V3)+  I]  I]  Am„am^{amxx)cos(KnXi) , 

m=l,3,5«=l,3,5  m=l,3,5  «=1,3,5 

00  00 

WU=-*2  X 

m=l,3,5  «=1,3,5 

00  00 

msin(«mx1  )cos(/r„x3 ) , 

m=l,3,5  n= 1,3,5 

OO  00  CO  CO 

W2>2  =  0.0598*2  Yj  Y  AmnalC°s(amX\  )COS(/r„*3)  +  0.0689x2  ^  £  A„A«C0S(«mX1)C0S(x-„X3), 

w=l,3,5  «=1,3,5  m=l,3,5  n=l,3,5 

00  00 

W2>3=-  X  Z^«^COS(^Xl)Sin(/t'«X3)) 

m=l,3,5  /2=1,3,5 


w3,i=-*2  X  X  A'""a"' v»sin(am*i  )sin( Kn  x3 ) , 

m=l,3,5  n=l,3,5 

00  00  00  00 

w3  2  =0.0860X2  X  X A»mamcos(amxi  )cos(*-„x3 ) - 0.267 1  x2  ^  ^Jm,I/£-);cos(«fflx1)cos(/£-/,x3) 

m=l,3,5  «=1,3,5  w=l,3,5  «=1,3,5 

00  00 

+  I  2>.,a  cos(«)),x1)sin(/t-(Ix3), 

m= 1,3,5  «=1,3,5 


^3,3=  *2  X  X'4'«''A'«C0S(r/'”A|  )C0S(AVfr)  ' 

w=l,3,5  n=l,3,5 

Once  the  displacement  gradients  are  determined,  they  will  be  substituted  into  the  perturbation 
integral  (14)  to  calculate  the  frequency  shift  caused  by  acceleration  [2].  In  our  program,  only 
second  and  third-order  elastic  constants  are  included  in  the  calculation,  all  other  nonlinear 
constants  such  as  third-order  dielectric,  piezoelectric,  electrostrictive  constants  are  ignored. 
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c KaLy  ~^KL^ay  +  cKaLNwy,N  +  cKNLywa,N 
+c KaLyAB^ AB  + ^AKaLyE A> 

-  _  0  /  c0  .  i  r?0 

eKLy  ~  eKLM  Wy,M  ~kKLyAB^AB  +bAKLyEA  ’ 
*KL  =bKLABSAB+%KLAEA ’ 


(14) 


bAKLy  =  bABCD  +  £0 SABSCD  ~  e0 (SACSBD  +  SAD$Bc)- 

Where,  <oM  is  the  unperturbed  angular  frequency,  «>  is  the  perturbed  angular  frequency, 
and  Acom  is  the  angular  frequency  shift.  Here,  cKaLy,  eKLy,  and  eKL  are  the  effective 


elastic,  piezoelectric,  and  dielectric  constants  respectively;  uy  is  a  specific  mode  shape 
function  in  the  unperturbed  condition  and  </>M  is  the  electrical  potential  for  this  specific  mode. 
,  ,  are  the  initial  stress,  strain  and  electrical  field  respectively  with  wr<N  defining 

the  displacement  gradient.  cKaly ,  dfAB,  sK,  arc  the  second-order  elastic,  piezoelectric  and 


dielectric  constants.  cKaLyAB,  kjKaLy ,  and  xkl/  are  the  third-order  elastic,  piezoelectric,  and 


dielectric  constants  respectively  and  bfALy  are  the  electrostrictive  constants.  s0  is  the 


permittivity  of  free  space. 
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Instruction 


Figure  1 .  A  Contoured  resonator  and  coordinate  system. 

Purpose: 

Calculating  free  vibration  frequencies  and  modes  of  AT-  or  SC-cut  bi-convex  or  plano-convex 
quartz  resonators,  either  unelectroded  or  fully  electroded. 

Input: 

2ho  —  resonator  center  thickness. 

R  —  electrode/resonator  mass  ratio=//2/z7(/;/?o).  p-electrode  mass  density.  2 h':  electrode 
thickness. 

Rc  —  radius  of  curvature  of  the  convex  surface  of  a  plano-convex  resonator,  or  half  of  the  radius 
of  curvature  of  the  convex  surface  of  a  bi-convex  resonator. 

n=  1,3,5,...  —  number  of  nodal  points  of  the  mode  along  the  resonator  thickness  direction  xi. 
m= 0,  2,  4  —  mode  order  along  the  in-plane  direction  xi  determined  by  the  Hermite  polynomial 
solutions  Hm(x\  )  of  the  modes. 

p= 0,  2,  4  —  mode  order  along  the  in-plane  direction  X3  detennined  the  Hennite  polynomial 
solutions  Hp(x 3)  of  the  modes. 

D  —  diameter  of  a  circular  resonator. 

2 L  and  2  W —  length  along  xi  and  width  along  X3  of  a  rectangular  resonator. 

Output: 

Frequency  (Hz). 

Mode  shape:  figure. 

Mode  shape:  data  array  (500  data  points  along  D,  or  2 L,  or  2  W). 

Note: 

The  frequencies  and  modes  of  AT-cut  resonators  are  calculated  using  Eqs.  (3.10)  and  (3.13)  of 

[1] .  The  frequencies  and  modes  of  SC-cut  resonators  are  calculated  using  Eqs.  (124)  and  (126)  of 

[2] .  The  modes  and  frequencies  in  [1]  and  [2]  are  for  resonators  unbounded  in  the  X1-X3  plane. 
Since  the  modes  die  out  quickly  away  from  the  center,  the  output  of  the  modes  is  calculated 
within  a  finite  circular  or  rectangular  region  predetennined  by  the  user  of  this  computer  code. 
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Instructions 
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sA 
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direction 

Figure  1.  A  half  space  of  AT -cut  quartz  and  coordinate  system. 

Purpose:  Calculating  the  displacement  components  (m  and  M3)  of  Rayleigh  waves  propagating  along  X3. 

Input: 

/ —  wave  frequency  (MHz)  or  X  —  wavelength  (mm). 


Output: 

X  —  wavelength  (mm)  or / —  wave  frequency  (MHz). 

Mode  variation  along  xi  (figures). 

Mode  variation  along  xo  (data  array,  501  data  points  along  X2). 


Theoretical  Background: 

For  Rayleigh  waves  propagating  along  the  X3  axis,  the  governing  equations  are: 


C22M,  22  +  2c24M2  23  +  C44M2  33  +  ^24^3,22  +  (C23  +  £44)113,23  +  £34  M 


'34“3,33 


:  PU2 


C24M2  22  +  (^23  £44)^2,23  ^34^2,33  £44^3,22  ^^34^3,23  ^33^3.33  3 


Consider  the  following  solution: 


u2  -  A  exp(£2x2) exp  i(k3x3  -  cot ) 
u3-B  exp(£2x2)  exp  i(k3x3  -  cot) 

Substituting  from  (2)  into  (1)  and  eliminating  the  conmion  exponential  factor,  we  obtain: 

(c22k\  +  2 c24ik2k3  -  c44k 2  +  pco2)  A  +  [c24C22  +  (£23  +  c44)ik2k3  -  c34k3\B  =  0 

[c24k2  +  (c23  +  c44)ik2k3  -  c34kl\A  +  (c44k22  +  2 c34ik2k3  -  c44k3  +  pco2  )B  -  0 


For  nontrivial  solutions,  the  determinant  of  the  coefficient  matrix  must  vanish,  thus 

(c22k2  +  2c24ik2k3  -  c44k3  +  pco2)  [c24k2  +  (c23  +  c44)ik2k3  -  c34k32] 
[c24k2  +  (c23  +  c44)ik2k3  -  c34k ;  ]  (c44k;  +  lc34ik2k3  -  c44k2  +  pco2) 


(1) 

(2) 

(3) 

(4) 


(4)  is  a  quartic  equation  for  k2 .  There  are  four  complex  roots.  For  surface  waves  the  displacements  must 
decay  as  X2  approaches  minus  infinity.  Only  the  two  roots  with  positive  real  parts  denoted  by  k‘2"  and  k{2) 
are  taken.  The  amplitude  ratios  au)  are  defined  by: 

B(i)  =  a(0A(i)  (5) 


where 


(6) 


a.(o  _  c22  ( K'f  +  2c24il(%  -  c44k;  +  par  ,=  ]  2 

c24(^2°)2  +  (c23  +  cM)ik"%  - c34k; 

The  displacement  field  can  be  expressed  as: 

u2  =  [A{  exp^'x,)  +  exp(Ar22’x2)]exp/(Ar3X3  -  cot) 
u3  =  [ a1AI  exp(k2l)x2)  +  a2A2  exp(k22>x2 )] exp i(k3x3  -  cot ) 

The  stress  components  needed  for  the  boundary  conditions  are  calculated  from: 

T22  —  c22u 2  t  +  c23u3  3  +  c24(ti2  3  +  u3  2) 

^23  —  ^24^2,2  ^34^3,3  ^44  (**2,3  **3,2  ) 

The  traction-free  boundary  conditions  are: 

T22  =  0 ,  r2,  =  0 ,  atx2  =0 

The  substitution  from  (7)  into  (9)  through  (8)  gives 

[(c22  +  c2Aal)k2)  +  (c24  +  Cyff,)/!:,]/!,  +  [(c22  +  c24a2)£22)  +  (c24  +  c23«2)//c3]d2  =  0 
[(c24  +  c44at)k2}  +  (c44  +  c34«1)il3]4  +  [(c24  +  c44a2)kf)  +  (c44  +  c34a2  )zl3  ]J2  =  0 

Nontrivial  solution  of  (9)  exists  if 

(c22  +  c24a[)^2n  +  (c24  +  c23a1  )ik3  ( c22  +  c24a2)k\2]  +  (c24  +  c23a2)ik3  ^ 

(c24  +  c^a^k™  +  (c44  +  c34at)ik3  (c24  +  c44a2)k2]  +  (c44  +  c34a2)ik3 


(7) 

(8) 

(9) 

(10) 


(11) 


(11)  is  an  equation  for  the  wave  speed  v  =  cot  k3 .  It  has  only  one  root.  The  corresponding  ratio  between 
A  i  and  A  2  is  determined  from  (10).  A\=l  is  used  in  the  output. 


